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I.  INTRODUCTION 

This  volume  provides  details  concerning  the  digital  computer 
code,  VANISH,  which  was  developed  for  the  Vulnerability  Analysis 
of  Nonlinear  Interacting  Soil  and  Hardened  Structures.  In  par- 
ticular, the  code  as  been  designed  to  be  efficient  and  reliable 
in  determining  the  response  of  a soil-covered  circular  arch  used 
as  an  aircraft  shelter.  However,  it  should  also  be  useful  in  ana- 
lyzing a variety  of  other  planar  structures.  Since  the  arch  of 
immediate  interest  is  composed  of  precast  reinforced  concrete  ribs, 
and  the  response  of  the  system  is  of  interest  all  the  way  to 
collapse,  special  measures  have  been  taken  to  account  for  nonlinear 
behavior  of  concrete  members,  and  of  nonlinear  behavior  of  the  soil. 

Volume  I of  this  report,  entitled  "Analysis  and  Experiment" 
gives  the  fundamental  concepts  underlying  the  VANISH  code.  In  par- 
ticular, Chapter  II  of  Volume  I gives  description  of  the  modeling 
techniques  to  be  used  for  various  structures,  derivations  of  the 
stiffness  and  mass  properties  of  the  three  types  of  finite  elements 
(beam,  bar  and  constant  strain  triangle  elements)  used  in  the  code, 
explains  how  the  code  handles  camping  and  applied  forces  (including 
blest  pressures),  and  presents  the  methods  used  for  incremental  static 
loading  and  for  dynamic  loading.  The  reader  will  find  understanding 
of  and  ready  reference  to  Chapter  II  of  Volume  I essential  to  under- 
standing of  the  latter  portions  of  the  present  volume,  which  document 
the  operation  of  the  code. 

Nevertheless,  most  of  the  knowledge  needed  merely  to  use  the 
VANISH  code  is  available  from  the  present  volume,  and  requires  only 
a general  understanding  of  how  the  code  actually  solves  a given 
problem. 

The  essential  steps  in  preparing  to  use  the  code  to  analyze  a 
structure  are:  1)  to  idealize  the  actual  structure  by  means  of 
straight  beam  elements  and/or  triangular  soil  elements,  2)  to  choose 
lines  of  possible  soil  cracking  between  triangular  elements. 
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if  such  nonlinear  behavior  of  the  soil  is  to  be  considered,  and 
introduce  short  bar  elements  across  these  potential  crack  lines, 

3)  to  number  all  of  the  nodes  and  elements  (by  element  type),  4)  to 
evaluate  the  stiffness  and  mass  properties  of  each  element,  5) 
to  choose  the  damping  coefficients  to  be  specified  for  one  or  two 
nodes,  6)  to  decide  what  loads  are  to  be  applied  (these  must  be 
concentrated  at  the  nodes),  and  7)  to  decide  what  results  are 
to  be  printed  out.  The  steps  in  describing  this  data  to  the 
computer  code  by  means  of  punched  cards  are  presented  in  Chapter 
II  of  this  volume. 

The  remaining  chapters  in  this  volume  give  detailed  infor- 
mation about  the  computer  code  itself,  so  that  potential  problems 
in  its  operation  and  extensions  of  its  capabilities  may  be  handled 
by  future  U'-'rs.  Chapter  III  gives  the  basic  function  and  Chapter 
IV  presents  a fr  'chart  of  each  subroutine  and  shows  how  all  of 
the  subroutines  are  related  to  one  another  (mostly  through  the 
MAIN  program).  Chapter  . ’ defines  the  important  variable  names 
used  and  identifies  the  s.oroutines  and  COMMON  blocks  where  they 
occur. 
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II.  INPUT/OUTPUT  INSTRUCTIONS 


A.  General 

The  first  step  in  the  analysis  of  a structure  by  the  finite 
element  method  is  to  divide  the  structure  into  finite  elements 
connected  at  points  called  nodes.  The  use  of  this  code  is  limited 
to  plane  stress  problems  employing  straight  beam  and  bar  elements 
and  plane  triangular  elements.  A suitable  x-y  coordinate  system 
should  be  selected  and  the  coordinates  of  the  nodes  determined. 

The  coordinates  of  the  nodes,  material  properties  of  the  elements, 
boundary  conditions,  damping  characteristics  and  loads  define  the 
problem  to  be  solved.  These  data  are  to  be  input  by  means  of  punched 
cards. 

B . Input  Cards  to  Define  the  Structure 

1 . Title  Card 

The  first  card  of  the  data  set  for  any  individual  problem 
must  have  a problem  number  in  columns  1 to  5.  This  is  to  enable 
a number  of  problems  to  be  solved  one  after  the  other.  A blank 
card  must  be  used  at  the  end  of  the  complete  data  set  to  specify 
the  end  of  all  the  problems.  The  title  card  also  allows  a title 
to  be  given  to  each  individual  problem  treated,  and  this  title  is 
printed  as  a heading  on  the  output. 

FORMAT  (15,  1 5A4) 

Columns  Data 

1 - 5 Problem  number. 

6-65  Title  of  the  problem. 

2 . Standard  Options  Card 

Standard  options  are  required  in  all  problems,  whereas  non- 
linear options  (following  card)  are  required  only  for  nonlinear 
problems. 
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FORMAT  (615) 


Co  1 umn  s 
1 - 5 

6 - 10 

11  - 15 


16  - 20 


21  - 25 


Data 

Parameter  specifying  the  type  of  problem. 

If  1,  static  analysis.  If  0 (or  blank), 
dynamic  analysis. 

Parameter  specifying  linear  or  nonlinear 
analysis.  If  1,  linear  analysis.  If  0 
(or  blank),  nonlinear  analysis. 

Parameter  specifying  whether  the  stresses 
are  to  be  determined.  If  1,  the  displace- 
ments only  are  determined.  If  0 (or  blank), 
the  displacements  and  stresses  are  determined. 

Parameter  controlling  printout  of  results. 

If  1,  the  displacements  (and  stresses  if 
calculated)  at  specified  nodes  only  are  printed. 
If  0 (or  blank),  displacements,  velocities, 
accelerations  (and  stresses  if  calculated) 
are  printed  for  all  the  nodes. 

Parameter  specifying  whether  dynamic  forces  are 
read  in  or  generated  internally  by  subroutine 
FORGEN.  If  1,  FORGEN  is  used.  If  0 or  blank 
forces  are  read  in. 


26  - 30  Parameter  specifying  inclusion  of  dead  loads. 

If  1,  dead  loads  area  neglected.  If  0,  dead 
loads  are  included. 


3.  Nonlinear  Options  Card 

This  card  is  omitted  if  the  analysis  is  linear. 


FORMAT  (615) 


Columns 
1 - 5 


6 - 10 


11  - 15 


Data 

Parameter  specifying  which  elements  (beams, 
bars,  or  both)  are  nonlinear.  If  0,  only 
bars  are  nonlinear.  If  1,  only  beams  are 
nonlinear.  If  2,  both  are  nonlinear. 

Parameter  specifying  type  of  beam  material 
in  axial  deformation.  If  0,  general  homo- 
geneous material  with  symmetric  yield 
levels.  If  1,  reinforced  concrete  (separate 
steel  and  concrete). 

Parameter  specifying  how  incremental  static 
loads  are  prescribed.  If  0,  increments  are 
set  internally.  If  1,  increments  are  read 
in . 
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16  - 20  Parameter  specifying  print-out  at  incremental 

static  load  levels.  If  0,  no  print-out  at 
each  load.  If  1,  standard  output  at  each  load. 

4.  P-'-'blem  Size  Card 


FORMAT  (515) 


Col umns 

1 - 5 
6 - 10 
11  - 15 
16  - 20 
21  - 25 


Data 

Number  of  nodal  points  (150  maximum*) 

Number  of  beam  elements  (30  maximum*) 

Number  of  bar  elements  (100  maximum*) 

Number  of  triangular  elements  (200  maximum*) 
Number  of  nodal  points  having  specified 
boundary  conditions  (33  maximum*) 

(*for  dimensions  in  present  code  version) 


5.  Nodal  Point  Coordinate  Cards 


One  card  is  required  for  each  nodal  point  in  numerical  sequence. 
FORMAT  (15,  5X,  2F10.0) 


Columns  Data 

1 - 5 Nodal  po'nt  number. 

6 - 10  Blank 

11-20  X-Coordinate  (units-feet) 

21  - 30  Y-Coordinate  (units-feet) 

6.  Beam  Element  Data  Cards  - Linear  Behavior 


These  cards  are  omitted  if  there  are  no  beam  elements  in  the 
finite  element  system.  One  c?"J  is  required  for  each  beam  element 
in  numerical  sequence. 

FORMAT  (315,  5X,  4F10.0) 


Columns  Data 

1 - 5 Beam  element  number. 

6-10  Nodal  point  number  of  one  end  of  beam  element. 

11-15  Nodal  point  number  of  the  other  end  of  the 

beam  element. 

16  - 20  Blank  2 

21  - 30  Area  of  cross-section  (in  ) ...  4. 

31  - 40  Moment  of  inertial  of  the  cross-section  (in  ) 

41  - 50  Modulus  of  elasticity  of  the  beam  material  (psi) 

51  - 60  Unit  weight  of  the  beam  material  (lb/ft  ) 
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I 
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7.  Beam  Element  Data  Cards  - Nonlinear  Behavior 


These  cards  are  omitted  if  the  beam  elements  remain  linear, 
but  are  needed  in  addition  to  the  linear  data  cards  for  nonlinear 
beams. 

a)  Nonlinear  Axial  Deformation  Properties 

(One  card  is  required  for  each  beam  element  in  numerical 
sequence.  Columns  41-80  may  be  left  blank  for  homogeneous 
beams ) . 


FORMAT  (15,  5X, 

7F10.0) 

Columns 

Data 

s 

1 - 5 

Beam  element  number 

6 - 10 

Blank 

11-20 

Yield  stress  (psi)  of  homogeneous  beam  or 
steel  of  R/C  beam 

21  - 30 

Rupture  strain  (in/in)  of  homogeneous  beam 
or  steel  of  R/C  beam 

Si 

31  - 40 

Rupture  stress  (psi)  of  homogeneous  beam  or 
steel  of  R/C  beam? 

41  - 50 

Area  of  steel  (in^)  of  R/C  beam 

51  - 60 

Compressive  yield  stress  (psi)  of  concrete 
R/C  beam  (f'). 

61  - 70 

Rupture  strain  in  compression  of  concrete  of 

5* 

R/C  beam 

71  - 80 

Rupture  stress  (psi)  in  compression  of  con- 
crete of  R/C  beam 

b)  Nonlinear  Bending  Deformation  Properties 

t 

(One  card  is 
sequence). 

required  for  each  beam  element  in  numerical 

FORMAT  (15,  5X, 

7F10.0) 

Columns 

Data 

1 

1 - 5 

Beam  element  number 

6 - 10 

Clank 

11-20 

Positive  yield  moment  of  each  rotary  spring 
(ft  - kips) 

21  - 30 

Absolute  value  of  negative  yield  moment  of  each 
rotary  spring  (ft  - kips) 

t 

31  - 40 

Positive  rupture  rotation  of  each  rotary 
spring  (radians) 
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41  - 50  Absolute  volume  of  rupture  rotation  of  each 

rotary  spring  (radians) 

51  - 60  Positive  rupture  moment  of  each  rotary  spring 

(ft  - kips) 

8.  Bar  Element  Data  Cards 

These  cards  are  omitted  if  there  are  no  bar  elements  in  the 
finite  element  system.  One  card  is  required  for  each  bar  element 
in  numerical  sequence. 

FORMAT  (315,  5X,  3F10.0,  110,  F10.0) 


Columns  Data 


1 - 

5 

Bar  element  number 

6 - 

10 

Nodal  point  number  of  one  end  of  the  bar 
element 

11  - 

15 

Nodal  point  number  of  the  other  end  of  the 
bar  element 

16  - 

20 

Blank  2 

21  - 

30 

Area  of  cros section  (in  ) 

31  - 

40 

Modulus  of  elasticity  of  the  bar  material  (psi) 

41  - 

50 

Unit  weight  of  the  bar  material  (lb/ft  ) 

51  - 

60 

Parameter  specifying  how  yielded  (tensile) 

stiffness  is  prescribed.  If  =1,  read  in 
by  next  parameter.  If  = 0,  specify  internally 
as  10  . 


61  - 70 

Yielded  (tensile)  stiffness  as 
pres  si ve  stiffness,  EA/L. 

ratio  to  com- 

9.  Trianqular 

Element  Data  Cards 

These  cards  are 

omitted  if  the  number  of  triangular  elements 

in  the- problem  is  zero.  One  card  is  required  for  each  triangular 

element  in  numerical 

sequence. 

FORMAT  (415,  4F10.0) 

Col  umns 

Data 

1 - 5 

Triangular  element  number 

6-10 

♦Node  number  of  1st  node  of  the 

triangular 

element 

11  - 15 

*Node  number  of  the  2nd  node 

16  - 20 

♦Node  number  of  the  3rd  node 

♦The  three  node  numbers  of  each  triangular  element  must  be  in 


counterclockwise  sequence. 


21  - 30  Thickness  of  the  element  (ft) 

31  - 40  Modulus  of  elasticity  of  the  material  (psi) 

41  - 50  Poisson's  ratio  of  the  material  « 

51  - 60  Unit  weight  of  the  material  (Ib/ft  ) 

10.  Boundary  Condition  Cards 

One  card  is  required  for  each  nodal  point  where  displacement 
boundary  conditions  are  prescribed.  The  code  is  written  such  that 
all  prescribed  boundary  conditions  must  have  zero  values.  These 
cards  need  not  be  in  any  numerical  sequence.  However,  the  total 
number  of  cards  must  be  equal  to  the  number  specified  on  the  problem 
size  card. 


FORMAT  (4.15) 

Columns 
1 - 5 

6 - 10 

11  - 15 

16  - 20 


Data 

Nodal  point  number  at  which  displacement 
boundary  conditions  are  specified. 

1,  if  the  u displacement  (displacement  in 
the  x direction)  is  zero.  Otherwise  blank. 

1,  if  the  v displacement  (displacement  in 
the  y direction)  is  zero.  Otherwise  blank. 

1,  if  the  rotation  is  zero.  Otherwise  blank. 


1 1 . Dynamic  Analysis  Parameters  Card 
FORMAT  (3F10.0,  110) 


Col umns 
1 - 10 

11  - 20 


21  - 30 
31  - 40 


Data 

Coefficient  by  which  to  multiply  the  stiffness 
matrix  to  get  the  "relative"  part  of  the 
damping  matrix. 

Coefficient  by  which  to  multiply  the  mass 
matrix  to  get  the  "absolute"  part  of  the 
damping  matrix. 

Magnitude  of  the,  time  step  (seconds). 

Maximum  number  of  time  steps  to  which  the 
analysis  is  to  be  carried  out. 


C. 


Input  Cards  to  Control  the  Output 

The  following  cards  may  be  omitted  if  the  displacements,  velo- 


cities, accelerations  and  stresses  are  desired  to  be  printed  out  for 


all  nodal  points.  Columns  15  - 20  in  the  "options  card"  (following 
the  title  card)  must  be  blank  in  this  case.  If  the  output  infor- 
mation is  desired  only  at  certain  (special)  nodal  points,  these 
nodal  point  numbers  are  input  using  the  following  cards.  Column  20 
} in  the  "options  card"  must  contain  1 in  this  case. 

1 . Number  of  Special  Nodes  Card 

FORMAT  (15) 

i 

} Col umns  Data 

1 - 5 Total  number  of  nodal  points  whose  displace- 

ments, stresses,  etc.  are  desired  to  be 
printed  out. 

t 2.  Cards  Specifying  Special  Node  Numbers 

FORMAT 

Columns  Data 

, 1 - 5 

6-10  Noda1  point  numbers  whose  displacements, 

- velocities,  etc.,  are  to  be  printed  out. 

61  - 65 
66  - 70 

& 

' Each  card  contains  14  numbers.  Use  as  many  cards  as  are  necessary 

to  specify  all  the  node  numbers. 

3.  Numbers  of  Special  Elements  Card 

* FORMAT  (515) 

Columns  Data 

Total  number  of  beam  elements  for  which 
member  forces  are  to  be  printed  out. 

Total  number  of  bar  elements  for  which 
for  which  forces  are  to  be  printed  out. 

Total  number  of  triangular  elements  for 
which  stresses  are  to  be  printed  out. 


1 - 5 
6 - 10 
11  - 15 
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4.  Cards  Specifying  Special  Element  Numbers 
FORMAT  (1415) 


Columns 

1 - 5 
6 - 10 


61  - 65 
66  - 70 


Data 


Element  numbers  whose  forces  (beams  and/or 
bars)  or  stresses  (triangles)  are  to  be 
printed  out;  all  beams  first,  then  all  bars, 
then  all  triangles. 
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D.  Input  Cards  for  Applied  Loads 

One  set  of  cards  must  be  included  for  each  time  step.  In 
each  set,  one  card  is  required  for  each  nodal  point  where  external 
loads  are  applied.  The  cards  need  not  be  in  any  particular  numerical 
sequence.  However,  the  last  card  in  a set  must  be  blank  indicating 
the  end  of  a set  of  loads  for  a particular  time  step. 

FORMAT  (15,  5X,  3F10.0) 

Columns  Data 

1 - 5 Node  number  where  the  external  loads  are 

applied. 

6 - 10  Blank 

11-20  X-Component  of  the  Load  (kips),  positive  upward 

21  - 30  Y-Component  of  the  load  (kips),  positive  to  right 

31  - 40  Moment  (ft  kips),  positive  counterclockwise 

Note  that  in  the  present  version  of  the  code,  there  is  no  pro- 

vision to  input  distributed  loads  (other  than  the  dead  weight  which 
is  automatically  calculated  using  the  unit  weights  of  the  materials). 

All  the  preceding  cards  specify  one  problem.  Data  for  a number 
of  problems  may  follow  one  after  the  other.  The  end  of  the  data  for 
all  problems  is  indicated  by  introducing  a blank  card  after  the 
complete  data. 
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E.  Output 

There  is  only  one  type  of  output,  namely  printed  output.  The 
following  information  is  printed  by  the  program  for  all  problems: 

a.  Input  data.  All  the  input  data  except  the  loads  are  printed 
before  the  execution  of  the  problem. 

b.  Error  messages.  The  code  checks  the  dimensions  provided 
against  the  size  of  the  problem  and  prints  out  error 
messages  if  the  size  exceeds  the  dimensions  provided. 

The  following  information  is  printed  for  each  load  increment, 

if  desired,  and  for  the  total  load  for  static  loading: 

a.  Displacements  of  each  node  or  designated  nodes. 

b.  Axial  force,  shear,  and  moment  at  each  end  of  each  beam 
element  or  designated  elements. 

c.  Axial  force  in  each  bar  element,  or  designated  elements. 

d.  Averaged  stresses  at  the  nodal  points  of  each  triangular 
element,  or  designated  elements. 

The  following  information  is  printed  for  each  time  step: 

a.  Nodal  displacements  and  velocities  of  each  node,  or  desig- 
nated nodes. 

b.  Forces  or  stresses  in  each  element,  or  in  designated 
elements,  if  called  for  (if  NOUT  =1). 
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III.  MAIN  PROGRAM  AND  SUBROUTINES 


A.  Main  Program 

The  MAIN  program  is  a control  program.  As  shown  in  the  accom- 
panying flowchart,  its  major  steps  are: 

1.  Read  and  print  the  problem  number,  title  and  options 
(static  or  dynamic  analysis,  linear  or  nonlinear  analysis'. 

2.  Call  subroutine  SETUP  in  which  most  of  the  data  are  read. 

3.  Call  subroutine  ASMBLE  which  assembles  the  stiffness  and 

mass  matrices. 

4.  If  the  analysis  is  static,  call  BCOND  to  incorporate  the 

boundary  conditions  into  the  stiffness  matrix,  call  DECOMP  to  decom- 

pose the  stiffness  matrix,  add  dead  and  applied  loads  if  appropriate, 
and  call  STSOL  to  solve  for  full  static  displacements. 

5.  If  the  analysis  is  static  but  may  be  nonlinear,  call  RESET 
to  check  for  yielding,  and  call  STNON  to  perform  an  incremental 
solution,  if  there  is  yielding. 

6.  If  the  analysis  is  dynamic,  perform  a static  solution  first 
if  dead  loads  are  considered,  then  perform  a dynamic  analysis,  using 
SOLVE  for  the  incremental  response  and  RESET  to  check  for  yielding 
and  modify  the  stiffness  matrix  and  restoring  force  vector  at  each 
time  step.  Dynamic  loads  may  either  be  read  in  or  calculated  by  FORGEN. 

The  subroutines  called  by  MAIN  are: 

SETUP 

ASMBLE 

BCOND 

DECOMP 

STSOL 

RESET 

STOON 

FORGEN 

SOLVE 
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B . Subroutine  SETUP 

This  is  the  primary  subroutine  for  reading  in  the  data,  although 
some  data  are  read  in  MAIN.  The  major  steps  in  this  subroutine  are: 

1.  Read  and  print  out  number  of  nodes,  number  of  elements  and 
other  quantities  specifying  the  size  of  the  problem. 

2.  Read  and  print  out  coordinates  of  the  nodes. 

3.  Read  and  print  out  data  for  the  beam  elements. 

4.  Read  and  print  out  data  for  the  bar  elements. 

5.  Read  and  print  out  data  for  the  triangular  elements. 

6.  Read  and  print  out  specified  displacement  boundary  conditions. 

7.  Read  and  print  parameters  for  dynamic  analysis. 

8.  Read  node  numbers  and  elements  whose  responses  are  to  be 
output. 

9.  Read  and  print  loads  for  static  analysis  (or  the  initial 
loads  for  dynamic  analysis). 

10.  Calculate  size  of  matrices. 

11.  Check  dimensions. 

SETUP  does  not  call  any  other  subroutines. 

C.  Subroutine  ASMBLE 

This  subroutine  is  for  assembling  the  stiffness  and  mass  matrices. 
The  modified  matrices  for  the  dynamic  analysis  (based  on  the  recursive 
form  of  Newmark's  3 method)  are  also  developed.  The  loads  due  to 
the  dead  weight  are  calculated  in  this  subroutine.  The  major  steps 
involved  in  this  subroutine  are: 

1.  Initialize  the  stiffness  and  mass  matrices. 

2.  Assemble  the  stiffness  and  matrices  ^or  the  beam  elements 
by  calling  subroutine  BEAM2. 

3.  Assemble  the  stiffness  and  mass  matrices  for  the  bar  elements 
by  calling  subroutine  BAR. 

4.  Assemble  the  stiffness  and  mass  matrices  for  the  triangular 
elements  by  calling  subroutine  CST. 
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5.  Store  the  stiffness  matrix  on  tapes  1 and  3. 

6.  Calculate  the  loads  due  to  the  deadweight. 

The  subroutines  called  by  ASMBLE  are: 

BEAM2 

BAR 

CST 

D . Subroutine  BCOND 

This  subroutine  incorporates  the  boundary  conditions  into  the 
stiffness  matrix,  that  is,  it  sets  the  off-diagonal  terms  equal  to 
zero  in  the  row  and  column  of  each  boundary  condition  and  it  sets  the 
diagonal  term  equal  to  one.  It  does  not  call  any  other  subroutine. 

E.  Subroutine  DECOMP 

This  subroutine  reduces  a symmetric  banded  matrix  into  an  upper 
triangular  matrix  by  Gauss  elimination.  There  is  no  provision  for 
pivoting  in  this  subroutine  since  the  matrices  involved  in  structural 
analysis  are  generally  positive  definite  and  diagonally  dominant.  No 
other  subroutines  are  called. 

F.  Subroutine  STSOL 

This  subroutine  is  used  in  the  static  analysis  only.  The  solution 
for  the  displacement  vector  is  obtained  using  this  subroutine,  which 
calls  subroutine  BSUB  for  this  purpose.  The  steps  involved  are: 

1.  Print  the  force  vector. 

2.  Incorporate  the  displacement  boundary  conditions  in  the  force 
vector. 

3.  Convert  the  forces  and  moments  to  inch  and  pound  units. 

4.  Call  subroutine  BSUB  to  solve  for  the  displacements. 

5.  Print  the  displacements. 

This  subroutine  calls  only  the  subroutine  BSUB. 
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G.  Subroutine  RESET 

This  subroutine  calculates  the  forces  and  stresses  in  individual 
structural  elements,  prints  these  data  if  desired,  checks  for  changes 
in  stiffness  or  restoring  force  due  to  yielding,  unyielding  or  rupture, 
and  modifies  the  stiffness  matrix  and  the  restoring  force  vector  as 
needed  for  the  next  time  or  load  increment.  The  major  steps  are: 

1.  Calculate  each  bar  force,  print  if  desired,  calculate  each 
bar  stiffness  for  next  step,  and  change  qlobal  stiffness  if  necessary. 

2.  Calculate  each  beam's  member  deformations  and  forces,  print 
forces  if  desired,  call  MODBM  to  check  for  stiffness,  energy,  or  restoring 
forces  changes  and  change  global  stiffness  if  necessary. 

3.  Calculate  each  triangle's  average  stresses  and  print,  if  neces- 
sary. 

Reset  calls  the  following  subroutines: 

BAR 

BEAMZ 

MODBM 

MULT2 

H.  Subroutine  STNON 

This  subroutine  sets  up  and  controls  the  incremental  static  solu- 
tion needed  when  thjre  is  nonlinear  behavior  under  static  loading.  For 
each  step  it  calculates  the  load  increment,  calls  STSOL  to  find  the 
increment  deflections,  adds  the  incremental  loads  and  deflections  to 
the  running  totals,  prints  the  results  if  desired,  and  calls  RESET  to 
check  for  (further)  nonlinear  behavior. 

STNON  calls  the  following  subroutines: 

BCOND 

DECOMP 

STSOL 

RESET 
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I.  Subroutine  FORGEN 


This  subroutine  allows  internal  calculation  of  nodal  forces  for 
dynamic  problems,  rather  than  having  to  read  them  in  at  each  time  step. 
It  is  especially  designed  to  handle  the  complex  loading  history  asso- 
ciated with  an  internal  blast  within  a soil-covered  circular  arch,  as 
simplified  from  scale  model  experiments.  It  does  not  call  any  other 
subroutines.  This  subroutine  was  developed  by  Mr.  Nash/DLYV/AFATL  at 
Eglin  AFB,  Florida  and  its  documentation  is  therefore  beyond  the  scope 
of  this  report. 

J.  Subroutine  SOLVE 

SOLVE  performs  a step-by-step  numerical  integration  of  the  dif- 
ferential equations  of  dynamic  equilibrium,  using  Newmark's  6 method 
with  8 = 1/4.  The  steps  involved  are: 

1.  Determine  the  accelerations  at  the  start  of  integration  from 
the  initial  conditions.  This  step  is  carried  out  for  the  1st  time 
step  only. 

2.  Add  the  dead  weight  to  the  force  vector  due  to  external  loads 
at  the  end  of  the  time  step.  (The  force  vector  is  input  in  the  MAIN 
program  for  each  time  step). 

3.  Derive  the  modified  force  vector  for  substitution  into  the 
recur; ive  equation  of  Newmark's  8 method, 

4.  Incorporate  the  displacement  boundary  conditions  into  the  modi- 
fied force  vector. 

5.  Solve  for  the  accelerations  at  the  end  of  the  time  step  by 
calling  subroutine  BSUB. 

6.  Determine  the  velocities  and  displacements  at  the  end  of  the 
time  step. 

7.  Store  the  above  values  as  starting  values  for  the  next  time 

step. 

8.  Print  the  desired  results. 

The  subroutine  SOLVE  calls  the  following  subroutines: 

BSUB 

MULT 
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K.  Subroutine  BEAMZ 

This  subroutine  derives  both  the  local  and  global  stiffness 
matrices  for  a straight  prismatic  beam  element  with  rotary  springs 
at  the  two  ends,  based  on  the  geometry  of  the  element  and  its  current 
stiffness  properties.  It  also  derives  the  compatibility  matrix  for  the 
beam  element  and  uses  it  to  convert  the  local  stiffness  evaluated  by 
subroutine  BMK  to  a global  stiffness.  The  subroutines  called  are: 

BMK 

MULTZ 

L.  Subroutine  BMK 

This  subroutine  derives  the  current  local  stiffness  matrix  for  a 
beam  element.  The  steps  are: 

1.  Determine  the  axial  stiffness  for  a general  beam  or  for  the 
steel  in  a reinforced  concrete  beam,  EA/L. 

2.  Add  in  the  concrete  axial  stiffness  if  the  the  beam  is  of 
reinforced  concrete. 

3.  Determine  the  current  bending  stiffness  based  on  the  current 
end  spring  stiffnesses. 

M,o  subroutines  are  called  by  BMK. 

M.  Subroutine  BAR 

This  subroutine  derives  the  stiffness  and  mass  matrices  of  a 
straight  bar  element  (axial  deformation  only).  The  stiffness  factors 
only  change  by  proportion  and  do  not  have  to  be  re-calculated  by  this 
subroutine  during  nonlinear  behavior.  The  subroutine  calls  no  other 
subroutines . 

N.  Subroutine  CST 

This  subroutine  derives  the  stiffness  and  mass  matrices  of  a con- 
stant strain  plane  triangular  element.  Each  node  is  considered  to  have 
two  translational  and  one  rotational  degree  of  freedom,  but  the  stiff- 
nesses corresponding  to  rotation  are  zero.  No  other  subroutines  are 
called. 
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O . Subroutine  BS'JB 

This  subroutine  solves  a set  of  simultaneous  linear  equations 
using  the  upper  triangular  matrix  obtained  from  subroutine  DECOMP. 

The  steps  involved  are: 

1.  Perform  a forward  substitution  for  the  input  vector  (forces). 

2.  Perform  a back-substitution  for  the  solution  vector  (dis- 
placements or  accelerations). 

No  other  subroutines  are  called. 

P.  Subroutine  MODBM 

This  subroutine,  which  is  called  from  RESET,  uses  HYST  to  cal- 
culate the  yielding  behavior  of  a single  beam  element  at  a particular 
step  of  the  solution.  It  considers  the  axial  deformation,  including 
both  steel  and  concrete  components  for  a reinforced  concrete  beam,  and 
the  bending  deformation  of  the  beam  at  each  end,  and  evaluates  the 
dissipated  energy  and  all  stiffness  and  restoring  force  changes.  It 
calls  the  following  subroutines: 

HYST 

BMK 

MULTZ 

Q.  Subroutine  HYST 

This  subroutine  determines  the  changes  in  a general,  nonsymmetric 
bilinear  hysteretic  diagram,  given  the  current  displacement,  permanent 
set  point  and  stiffness.  The  changes  include  possible  stiffness 
changes  due  to  yielding  or  unyielding,  possible  restoring  force  drop  due 
to  rupture,  and  energy  changes  due  to  any  (further)  yielding.  The 
major  steps  are: 

1.  Determine  if  positive  or  negative  rupture  has  occurred 
and  modify  the  restoring  force  and  stiffness  accordingly. 

2.  Determine  if  the  system  remains  elastic  or  yielded  in  either 
direction  or  yields  for  the  first  time,  and  set  the  stiffness  keys 
accordingly. 
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3.  Calculate  the  dissipated  energy  at  the  current  displacement 
as  an  incremental  part  associated  with  the  permanent  set. 

No  other  subroutines  are  called  by  HYST. 

R.  Subroutine  MULT 

This  subroutine  performs  a matrix  vector  multiplication  of  the 
form  y = A x where  A is  a symmetric  and  banded  matrix  one 
half  of  which  is  stored  as  a singly  subscripted  array.  No  other 
subroutines  are  called. 

S . Subroutine  MULT2 

This  subroutine  performs  a matrix  multiplication  of  the  form 
[A]  = [B]  [C],  where  [B]  has  dimension  N x N1 , [C]  dimensions 
N1  x M,  and  [A]  comes  out  with  dimensions  N x M.  No  other  sub- 
routines are  called. 
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IV.  FLOWCHARTS 

For  the  MAIN  program  and  for  each  of  the  subroutines  described 
briefly  in  Chapter  III,  a flowchart  giving  more  details  is  presented 
in  the  following  pages.  Some  of  the  variables  used  are  defined  in 
the  flowcharts;  others  may  be  found  in  the  descriptions  of  major 
variables  to  be  found  in  Chapter  VI. 
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A.  MAIN 


START 


Read  Problem  Number,  Title 


Prob.  No.  = 0?. 


Print  Problem  Number,  Title 


Read  and  Print  Standard  Options: 
NTYPE,  NUN,  NRES,  NOVT,  NFORC , NDEAD 


Nonlinear  Problem! 


Read  and  Echo  Print  Nonlinear  Options 
NLEXT,  NLAX,  ISLOAD,  ISPRT,  NDEAD 


Call  SETUP  (to  read  remainder  of 
input  data  except  dynamic  forces) 


Call  ASMBLE  (to  get  linear  static 
stiffness  matrix,  Ak,  on  tape  1) 


Initialize  to  zero  for  this  problem: 
NS,  NT,  displacement,  velocity, 
acceleration  and  force  vectors 


A.  MAIN  (Cont'd.) 


Dynamic  problem 
.rieglecting  dead  loadsjt- 


Read  AK  from  Tape  1 


Call  BCOND  (to  incorporate 
boundary  conditions  into  AK) 


Call  DECOMP  (to  decompose 
AK  for  equation  solving) 


Select  full  static  load,  FF, 
(with  or  w/o  dead  loads)  and 
print  headings  for  static  solution 


Call  STSOL  (to  obtain  and  print  ] 
linear  solution  at  full  static  load)  | 


Call  RESET  (to  check  for  nonlinear 
response  and  indicate  with  KCHAK  ii 
an  incremental  static  solution 
is  needed) 


KCHAK  = 0?„ 


NT  * NT  + 1 


Call  STNON  (to  perform 
incremental  static  solution) 
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Read  and  print  load  increments 
for  static  nonlinear  analysis 

E 

Calculate  half bandwidth  \ 
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G.  Subroutine  RESET 
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friy  beam  elements?" 


(For  each  beam  element) 


Form  current  element  stiffness  and 
compatibility  matrices,  using  subroutine  BEAM2 


Call  MULT2  to  get  incremental 
member  deformations,  DUBM=B*UF 


Call  MULT2  to  get  incremental 
member  forces,  DPBM  = $BM*DUMB 


Add  incremental  member  forces  to  totals, 
PBM=PBM  + DPBM,  and  convert 
to  inch  and  pound  units 


Print  beam  forces,  if  desired 


Nonlinear  Beams? 


Call  MULT2  to  get  current  total 
member  deformations,  UBM=B*UI 


Call  MODBM  to  check  for  yielding, 
unyielding  or  rupture,  to  find  the  stiffness, 
energy  dissipation  and  force,  and  to 
print  member  forces,  if  desired 


BB 


Initialize  KCHAK,  response  quantities 
(UI,  UF,  VF,  FR  and  PGM)  and  nonlinear  quantities 
(ENDIS,  KYBMS,  KYBMC,  LPS,  EPSC,  KYBM, 

HRPS,  HR,  and  KYBR)  to  zero 


Calculate  next  load  increment,  AF, 
from  percentage  of  FF 


Call  STSO 
deformations. 


Add  incremental  loads  and  delections 
to  running  totals,  VF  and  UI 


Call  Reset  to  check  for  yielding 
and  rupture  and  modify  AK,  if  necessary 


KCHAK=0? 


Read  AK  from  tape  3,  incorporate 
boundary  conditions  and  decompose 
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K.  Subroutine  BEAM2 


START 


Initialize  global  stiffness 
matrix,  S,  to  zero 


Calculate  geometric  parameters 
(length,  cos  and  sin  of  angle 
with  x axis,  a) 


Form  compatibility  matrix,  BB 


Form  transpose  of  BB,  BT 


Calculate  axial  factor  AE/L 
and  bending  factor  EI/L,  then 
call  BMK  to  get  3 x 3 local 
stiffness  matrix  SBM 


Call  MULT2  twice  to  get  stiffness 
matrix  in  global  coordinates, 
S=BT*SBM*BB 


RETURN 


L.  Subroutine  BMK 


START 


Set  second  and  third  terms 
of  first  row  and  column  of 
SBM  equal  to  zero 


Set  general  or  steel  portion  of  axial 
stiffness  equal  to  the  elastic 
yielded,  or  ruptured  value, 
according  to  the  key  KYA 


Concrete  beam? 


j Yes  • 

Set  concrete  portion  of  axial 
stiffness  equal  to  the  elastic, 
yielded  or  ruptured  value, 
according  to  the  key  KYAC 


for  each  end  of  the  beam 
1 


l 

Set  rotary  end  spring  stiffness  equal  j 
to  the  elastic,  yielded,  or  ruptured  i- 
value,  according  to  the  key  KYB1  j 


Form  the  remaining  terms  of  the  SBM 
matrix  (22,  23,  32  and  33  elements) 


P.  Subroutine  MODBM 


[OVal 


Initialize  dissipated  energy,  ENDI  and  load 
drop  due  to  rupture,  PDROP,  to  zero 


Calculate  remaining  basic  diagram  parameters 
SY,  FYP,  FYN,  FRP,  FRN,  CXYP,  CXYN 


Positive  rupture? 
^ (X  > XRP?1^ 


Negative  rupture? 
(X  < XRN?)  . 


PDROP =- FRN 


Increment  rupture  key, 
NEWR,  and  set  current 
stiffness  key,  KY,  to  Z 


Calculate  differences  between  current 
deformation  and  current  yield  points: 
DXP=X-CXYP 
DXN=X-CXYN 
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I.  Subroutine  CST. 


i 
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M.  Subroutine  STSOL. 


statements.  The  same  variable  names  are  used  in  the  MAIN  program 
and  the  major  subroutines  called  Dy  MAIN,  namely,  subroutines 
t SETUP,  ASMBLE,  STSOL,  RESET,  STNON,  and  SOLVE. 


Variables  in  the 

First  COMMON  Block  (Inteqer,  Dimensions  and 

Options ) 

Variable  Name 

Description 

NNOD 

Number  of  nodes 

NBMEL 

Number  of  beam  elements 

NBREL 

Number  of  bar  elements 

NTREL 

Number  of  triangular  elements 

NROW 

Number  of  rows  of  the  stiffness  matrix  and 
mass  matrix 

NHB 

Half  band  width  (including  diagonal  of  the 
stiffness  matrix 

NBOU 

Number  of  nodes  with  prescribed  displace- 
ment boundary  conditions 

IBOU 

Total  number  of  prescribed  boundary  conditions 

NSMAX 

Maximum  number  of  time  steps  for  dynamic 
analysis 

NMAX 

Available  computer  storage  storage  for  the 
stiffness  matrix.  (The  value  of  NMAX  depends 
upon  the  computer  and  the  dimensions  of  other 
variables.  This  value  should  be  specified 
in  the  MAIN  program. 

I ERR 

Error  flag 

I ERR 

NATURE  OF  ERROR 

101 

Number  of  beam  elements  exceeds  dimensions 
provided 

102 

Number  of  bar  elements  exceeds  dimensions 
provided 

103 

Number  of  triangular  elements  exceeds 
dimensions  provided 

s 

1 


ik  t-v.. ■ 


in 


104 

Number  of  rows  of  the  matrices  exceeds 
dimensions  provided 

105 

Matrix  size  is  greater  than  NMAX 
(available  core  storage) 

106 

Number  of  boundary  conditions  exceeds 
dimensions  provided 

NTYPE 

Option  parameter  to  specify  static  or  dy- 
namic analysis 

NUN 

Option  parameter  specifying  linear  or  non- 
linear analysis 

NRES 

Option  parameter  specifying  whether  or  not 
stresses  are  to  be  determined  at  each  step 

NOUT 

Option  parameter  controlling  the  output 

NNDO 

Total  number  of  special  nodes  where  responses 
are  chosen  to  be  printed  out 

NNBL 

Total  number  of  special  beam  elements  whose 
forces  are  to  be  printed  out 

NNBR 

Total  number  of  special  bar  elements  whose 
forces  are  to  be  printed  out 

NNTR 

Total  number  of  special  triangular  elements 
whose  stresses  are  to  be  printed  out 

NT 

Counter  for  nonlinear  static  load  increments 

B.  Variables  in  the  Second  COMMON  Block  (Dynamic  Parameters) 


Variable  Name 
Cl 

C2 

H 

BETA 


Description 

Constant  to  multiply  the  stiffness  matrix 
to  get  one  part  of  the  damping  matrix 

Constant  to  multiply  the  mass  matrix  to  get 
the  other  part  of  the  damping  matrix 

Magnitude  of  the  time  increment  for  dynamic 
analysis 

Newmark's  s parameter 


C.  Variables  in  the  Third  COMMON  Block  (Element  Node  Numbers) 


Variable  Name  Description 


NDBM  Beam  element  node  numbers 

NDBR  Bar  element  node  numbers 
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Triangular  element  node  numbers 
Special  nodes  for  output 
Special  beam  elements  for  print-out 
Special  bar  elements  for  print-out 
Special  triangular  elements  for  print-out 
Specified  boundary  conditions 


Variables  in  the  Fourth  COMMON  Block  (Vectors  of  Linear  System 
Properties 


Variable  Name 


Description 

x coordinates  of  nodes,  positive  to  right 

y coordinates  of  nodes,  positive  up 

Moduli  of  elasticity  of  beam  elements 

Cross-sectional  areas  of  beam  elements 

Moments  of  inertia  of  beam  elements 

Unit  weights  of  beam  elements 

Moduli  of  elasticity  of  bar  elements 

Cross-sectional  areas  of  bar  elements 

Unit  weights  of  bar  elements 

Moduli  of  elasticity  for  triangular  elements 

Poisson's  ratios  for  triangular  elements 

Thicknesses  of  triangular  elements 

Lin’-t  weights  of  triangular  elements 

Stiffness  matrix 

Mass  matrix 

Disolacements  at  the  beginning  of  a time  sfec 
(also  displacements  in  stat'c  analysis} 

Vr\  * q r ' i * r rj  + {"*0  0“  PH  ~ Q’f  3 !“"■  o vpr^ 

Accelerations  at  the  beginning  of  a time  step 

Forces  at  the  beginning  of  a time  step 
(also  forces  for  static  analysis) 

Displacements  at  the  end  of  a time  step 

Velocities  at  the  end  of  a time  step 

Accelerations  at  the  end  of  a time  step 
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FF  Forces  at  the  end  of  a time  step 

OW  Dead  weight  forces 

S Stiffness  matrix  of  an  element.  (Same 

variable  name  is  used  for  the  beam,  bar 
and  triangular  elements) 

SM  Mass  matrix  of  an  element  (Same  name  for 

all  elements) 

PBM  Forces  and  moments  in  beam  elements 

PBR  Axial  forces  in  bar  elements 

PTR  Stresses  in  triangular  elements 

PBBM  Beam  element  forces  and  moments  due  to 

dead  load 

PBBR  3ar  element  forces  due  to  dead  load 

PBTR  Triangle  element  stresses  due  to  dead  load 


f 


* 


I 


E.  Variables  in  the  Fifth  COMMON  Block  (Nonlinear  Potions  and 
Parameters ) 


Variable  Name  Description 


NLEXT 

NLAX 

ISPRT 

ES 

CRACK 

NPCTLD 

BL 

NDEAD 

EN 


Option  parameter  to  specify  if  beams,  bars 
or  both  are  nonlinear 

Option  parameter  to  specify  homogeneous 
or  reinforced  concrete  beams 

Option  parameter  allowing  reading  or 
interval  prescription  of  static  load  incre 
ments 

Modu7us  of  elasticity  of  steel 

Percentage  of  compressive  yield  load  at 
tensile  cracking  of  concrete 

Number  of  incremental  static  load  percentages 
soeci tied 

Beam  or  bar  length 

Option  parameter  specifying  whether  dead 
loads  are  included  or  neglected 

Total  dissipated  energy  in  a structure 
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F . Variables  in  the  Sixth  COMMON  Block  (Nonlinear  Element  Properties 

and  Keys -Vectors  and  Matrices! 

Variable  Name  Description 


EYS 

ERS 

SSBMS 

AS 

EYC 

ERC 

SSBMC 

HS 

HRYP 

HRYN 

HRRP 

HRRN 

SSBM 

EPS 

KYBMS 

EPSC 

KYBMC 

HRPS 

KYBM 

ENDIS 

SSBR 

FR 

HR 

KYBR 


Yield  strain  of  steel  in  a R/C  beam 
Rupture  strain  of  steel  in  a R/C  beam 
Second  slope  ratio  of  beam  steel 
Area  of  beam  steel 

Yield  strain  of  concrete  in  a R/C  beam 

Rupture  strain  of  concrete  in  a R/C  beam 

Second  slope  ratio  of  beam  concrete 

Rotary  hinge  stiffness 

Positive  yield  value  of  hinge  rotation 

Negative  yield  value  of  hinge  rotation 

Positive  rupture  value  of  hinge  notation 

Negative  rupture  value  of  hinge  rotation 

Second  slope  ratio  of  beam  hinge 

Strain  at  permanent  set  point  of  steel 

Yield  key  for  beam  steel 

Strain  at  permanent  set  point  of  concrete 

Yield  key  for  beam  concrete 

Hinge  rotation  at  permanent  set  point 

Yield  key  for  beam  hinge 

Dissipated  energy  in  a beam 

Second  slope  ratio  for  a bar 

Restoring  force 

Hinge  rotation 

Yield  key  for  a ban 
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6 . Variables  in  the  Seventh  COMMON  Block 
Variable  Name  Description 


PCTLD 

SBM 

BB 

BT 

VBM 

VBR 


Percentage  of  full  static  load  in  an 
incremental  static  solution 

Local  stiffness  matrix  of  a beam  element 
(3x3) 

Compatibility  matrix  of  a beam  element 
(3  x 6) 

Transpose  of  BB 

Local  deformations  of  a beam  element 
Local  deformation  of  a bar  element 


H.  Variables  not  included  in  COMMON  statements: 

Additional  variables  used  in  the  MAIN  program  and  the  subroutines 
are  described  below. 


MAIN 


Variable  Name 


Description 


TITLE  Problem  title 

NPROB  Problem  number 

NS  Time  step  number  for  dynamic  analysis 

2.  Subroutine  SETUP 


Variable  Name 


Description 


Form 
IBC,  BC 
3.  Subroutine  ASM3LE 


Variable  used  for  obiect-time-format  for  print- 
out of  prescribed  boundary  conditions 

Work  areas  used  for  reading  in  boundary  condi- 
tions and  loads. 


Variable  Name 
KK,  XX,  YY 


Description 

Work  areas  used  for  assembling  the  stiffness  and 
mass  matrices 
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4. Subroutine  BCOND 


Variable  Name 

Description 

A 

Matrix  to  have  boundary  conditions  incorporated 
(stored  as  a singly  subscripted  variable) 

NROW 

Number  of  rows  of  A 

M 

Half  band  with  (including  diagonal)  of  A 

NA 

NROW  * M = Total  storage  for  A 

5.  Subroutine 

DECOMP 

Variable  Name 

Description 

A 

Matrix  to  be  decomposed 

NROW,  m,  m 

As  in  Subroutine  BCOND 

6.  Subroutine 

STSOL 

Variable  Name 

Description 

A,  NROW,  M,  NA 

Same  as  in  subroutines  BCOND  and  DECOMP 

X 

Solution  vector  (Force  vector  on  entering 
the  subroutine,  Displacement  vector  when 
coming  out  of  the  subroutine). 

F 

Force  vector  corresponding  to  the  static 
loads 

7.  Subroutine 

RESET 

V ariable  Name 

Description 

SC 

Scale  factor  used  for  converting  from  inches 
and  pounds  to  feet  and  kips 

XX,  YY,  UUU,  FF  Work  areas  for  determining  the  stresses  in  the 

elements 

8.  Subroutine 

STNON 

Variable  Name 

Description 

IE 

Index  for  two  end  springs  of  a beam 

NAK 

Total  size  of  AK 

P 

Percentage  of  full  static  load 

9.  Subroutine 

FORGEN  (To  be  included  later  — in  final  report) 

10.  Subroutine  SOLVE 


Variable  Name 
HH,  DTI,  DT2 


Description 


Parameters  in  Newtr.ark  Beta  method 


11.  Subroutines  8EAM2 


Variable  Name 


Descri ption 


MS3,  MS6,  POSTS,  Temporary  work  SDaces  for  matrix  multiplications 
PRES  (BE 


SI,  CO 


(BE 

Sine  and  cosine 

x coordinates  of  the  nodes 

y coordinates  of  the  nodes 

Modulus  of  elasticity 

Area  of  cross  section 

Moment  of  inertia  (only  in  BEAM2) 

Unit  weight 


12.  Subroutine  SMK 
Variable  Name  Description 

KYA,  KYB,  R,C,  SBM  Local  names  for  arguments,  stiffness 
ASC,  KYAC,  SAC,  BS, 

KYB1 , KYB2,  SB, 


13.  Subroutine  BAR 
Variable  Name 


Description 


XX,  YY,  BE,  BA,  BW  Same  as  in  BEAM2 


14.  Subroutine  CST 
Variable  name 


Description 


x coordinates  of  the  nodes 
y coordinates  of  the  nodes 
Thickness  of  the  element 
Modulus  of  elasticity 
Poisson's  ratio 
Unit  weight 
Area  coordinates 


^ •; 


1 5.  Subroutine 
V ariable  name 
A,  NROW,  M,  NA 


Description 


Same  as  in  subroutine  DECOMP 


16.  Subroutine  MODBM 


Variable  Name 


AO,  AOC,  A1 
KCHFR 

17.  Subroutine 
Variable  Name 
DXN,  DXP 

ENDI , ENXPS 
CXYN,  CXYP 
FRP,  FRN 
KY 


PDROP 
S,  SS,  SY 


XRP,  XP.N 


XYP,  XYN 

17.  Subroutine 


Variable  Name 


Description 

Increment  in  hinge  rotation 
Current  axial  and  bending  stiffnesses 
Indicator  of  rupture 


Descriptions 

Differences  from  negative  and  positive  yield 
points 

Dissipated  and  permanent  set  energies 

Current  positive  and  negative  yield  points 

Positive  and  negative  rupture  forces 

Current  stiffness  key 

Counter  on  rupture 

Counter  on  yield 

Loss  of  force  due  to  rupture 

Stiffness  parameters 

Displacement 

Permanent  set  displacement 

Positive  and  negative  rupture 
displacments 

Positive  and  negative  yield  disolacements 


De^c*"-!  on  on 

Matrix 

Vector  to  be  multiplied  with  A 
Resultant  vector  (Y)  = [A]  (x) 

Number  of  rows  in  A 

Half  band  width  (including  diagonal)  of  A 
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